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Abstract 

Background: Resource availability and predation are believed to affect community dynamics and composition. 
Although the effects of resource availability and predation on prey communities are usually studied in isolation, 
these factors can also have interactive effects, especially since the outcome of competition under shared predation 
is expected to depend on resource availability. However, there are few experimental studies that test the interactive 
roles of resources and predation on dynamics of more complex multispecies communities. Here, we examine the 
importance of competition and predation on microbial community dynamics in a resource pulse environment. 

Results: We manipulated resource availability and predation simultaneously in a microbial microcosm experiment, 
where a bacterial community was exposed to the protozoan predator Tetrahymena therrnophila in three different 
resource concentrations (low, intermediate and high). The prey community consisted of three heterotrophic 
bacterial species: Bacillus cereus, Serratia marcescens and Novosphingobium capsulatum, all feeding on a shared plant 
detritus medium. In fresh culture media, all species grew in all resource concentrations used. However, during 
experiments without any addition of extra resources, the existing resources were soon depleted to very low levels, 
slowing growth of the three bacterial species. Prior to the microcosm experiment, we measured the competitive 
ability and grazing resistance, i.e. reduced vulnerability to predation, of each prey species. The three species differed 
in allocation patterns: in general, N. capsulatum had the best competitive abilities and B. cereus had good grazing 
resistance abilities. In the long-term microcosm experiment, N. capsulatum dominated the community without 
predation and, with predation, B. cereus was the dominant species in the intermediate and high resource 
environments. 

Conclusions: Short-term, single-species assays revealed significant differences in the allocation of competitive and 
defensive traits among the prey species. Based on these differences, we were, to some extent, able to predict how 
the long-term community structure, e.g. species dominance, is modified by the resource availability and predation 
interaction in pulsed resource environments. Our results are consistent with theoretical predictions and also 
highlight the importance of interactive effects of resource competition and predation, suggesting that these factors 
should not be studied in isolation. 
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Background 

The resource competition theory predicts that competi- 
tion will reduce species abundance and ultimately lead 
to dominance by a single species in cases where there is 
a constant supply of a single limiting resource [1]. How- 
ever, resource availability is seldom constant in natural 
environments; instead, the frequency and magnitude of 
resource supply often varies resulting in the alternation 
of low and high resource availability [2,3]. Perhaps the 
most drastic changes in resource availability are caused 
by resource pulses. Resource pulses are defined as low 
frequency temporal events leading to increased resource 
availability with a short duration and a large magnitude [4]. 
The temporal variability in resource supply during resource 
pulse events makes predicting the outcome of resource 
competition and community dynamics more challenging. 
One way to approach this problem is to investigate the 
characteristics related to resource use among the individual 
species in the community. In general, a species' ability to 
grow can differ between low and high resource environ- 
ments so that the competitive dominance between species 
changes as a function of the resource concentration. In en- 
vironments where resource pulses occur, rapid increases in 
resource concentrations can favor high resource specialists 
since their densities can respond rapidly to increasing re- 
sources. In contrast, long periods of low resource concen- 
tration can favor species that are able to maintain positive 
growth rates with very limited resources. Therefore, tem- 
poral changes in resource availability may allow different 
species to co-exist in the communities [1,5]. 

Predation, in addition to resource availability, is a major 
force shaping community structure. Predation often greatly 
reduces prey population density and alters community 
composition and species diversity [6,7]. Predation can have 
a positive effect on prey community diversity when preda- 
tors feed more on superior competitors, which would, 
without predators, dominate the community [8-10]. Fur- 
thermore, many studies have shown that grazing resistance 
is often costly and the cost is seen as reduced competitive 
ability [11,12]. The effect of resources and predation on 
prey community is often linked. The effect of predation on 
prey community composition depends on resource avail- 
ability [13-16]. Using a theoretical approach, Holt et al. 
(1994) found that traits related to resource competition are 
more important for prey fitness when resources are scarce. 
However, when resources are more abundant, traits related 
to grazing resistance can play a more significant role in 
prey fitness [16]. 

The aim of this study was to investigate the relative 
roles of resource competition and predation on prey com- 
munity dynamics after a resource pulse. In an earlier study 
[17] with the same system, fluctuating resource supply and 
predation allowed different prey species to prevail under 
different fluctuation regimes through trade-offs between 



competitive ability and grazing resistance. The present 
study continued this work by investigating the community 
dynamics in an environment where resources do not fluc- 
tuate but are provided as a single large-scale pulse and are 
thereafter depleted by the bacterial community. We tested 
how the amount of resources (intensity of the resource 
competition) and predation by the protozoan predator T. 
thermophila affect the population dynamics and compos- 
ition of an aquatic bacterial prey community. To answer 
these questions, we first conducted a short-term growth 
and feeding experiment to estimate the competitive ability 
and grazing resistance of each prey species. All measure- 
ments were conducted in three resource concentrations, 
low, intermediate and high. We used short-term popula- 
tion growth rate as an estimate of competitive ability and 
grazing resistance as an estimate of biomass reduction by 
predators. To investigate long-term community dynamics, 
we then conducted a 21-day long microcosm experiment 
where a three-species prey community was cultured in 
three resource concentrations with and without predators. 
During the experiment, we monitored population densities 
and prey community composition with direct microscope 
counts and a dilution plating method. The experimental 
set up simulated a resource pulse and we monitored the re- 
sultant community dynamics. We predicted that resources 
and predation should have a substantial effect on the prey 
community in the microcosm experiment, based on theory 
[6,7] and previous studies of similar microbial systems 
[9,10,17-19]. We also predicted that traits of individual spe- 
cies should explain community dynamics with traits related 
to competition being more important in the absence of 
predators and traits related to grazing resistance being 
more important when predators are present. In our case, 
we concluded, based on competitive ability and grazing re- 
sistance measurements, that N. capsulatum was the best 
resource competitor in low and intermediate resource con- 
centrations and B. cereus had overall good grazing resist- 
ance abilities. The role of predation is considered to be less 
important in low resource environments where resource 
competition is more important [16]. One reason behind 
this prediction is that in less productive environments 
predator population sizes are also smaller and thus preda- 
tion pressure is lower. The microcosm community experi- 
ment confirmed our predictions and the prey species that 
was a better resource competitor was dominant in the ab- 
sence of the predator and the more grazing resistant spe- 
cies dominated in the presence of predation. Furthermore, 
the role of predation was less significant in the low re- 
source environment where traits related to competitive 
ability seem to be more important in determining the iden- 
tity of the dominant species. 

Results 

For the simplified results summary see Table 1. 
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Table 1 Summary of the results from the trait and community experiments 

Predator absent Predator present 

Resource level Best competitive ability Dominant species in Best grazing resistance Dominant species in 

the community experiment the community experiment 

Low N.c N.c B.c/S.m N.c 

Intermediate N. c N. c B.c/N.c Be 

High S.m N.c B.c/N.c B. c 

Species codes: B. cereus (6. c), S. marcescens (S. m) and N. capsulatum (N. c). 



Trait measurements 
Competitive ability 

Species growth rates increased with increasing resource 
concentration, but the competitive rank between species 
varied among resource concentrations (Figure 1A; re- 
source treatment: F 2 .98 = 85.7; p < 0.0001, species iden- 
tity x resource treatment: F 6 . 98 = 36.7; p < 0.0001). We 
found that N. capsulatum had a higher growth rate than 
the two other species in low and intermediate resource 
concentrations when comparing species growth rates in 
different resource concentrations (Figure 1A; homogenous 
subsets; low resource concentration: a: p = 0.10; b: p = 
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Figure 1 Performance of the prey species in resource 
competition and predation environments. Fitness measured as 
population growth rate (A) and grazing resistance (B) (mean ± S.E.). 
In both panels letters (a, b or c) indicate subsets that are statistically 
not different from each other. 



1.00; intermediate resource concentration: a: p = 0.43; b: 
p = 1.00). S. marcescens had the highest growth rate, N. 
capsulatum had the second highest and B. cereus had 
the lowest, in high resource concentrations (Figure 1A; 
all species in different homogenous subsets). 

Grazing resistance 

Grazing resistance, as with population growth rates, was 
different among species and rank among species depended 
on the resource concentration (Figure IB; species identity: 
F3.98 = 79.4; p < 0.0001, resource treatment: F 2 . 98 = 19.5; 
p < 0.0001). The grazing resistance of S. marcescens de- 
creased with increasing resource concentration, the graz- 
ing resistance of N. capsulatum increased, and resource 
concentration had no effect on the grazing resistance of 
B. cereus (Figure IB, species identity x resource treatment: 
F 6 . 98 = 75.2; p < 0.0001). We found that N. capsulatum had 
a lower grazing resistance in low resource concentrations 
than the two other species (Figure IB; homogenous subsets; 
a: p = 0.84; b: p = 1.00). In intermediate and high resource 
concentrations, S. marcescens had lower grazing resistance 
than the other species (Figure IB, homogenous subsets; 
intermediate resource concentration; a: p = 0.13; b: p = 1.00; 
high resource concentration; a: p = 0.12; b: p = 1.00). 

Community experiment 
Population densities 

The total number of prey individuals was lower when 
predators were present in all resource concentrations and 
each concentration produced low, intermediate and high 
population densities (Figure 2; predator treatment: Fi i8 = 
413; p < 0.0001, resource treatment: F 218 = 440; p < 0.0001, 
all in different homogenous subsets). Predator population 
density was lower in low resource concentrations com- 
pared to intermediate and high resource concentrations 
(Figure 2; resource treatment: F 2 .9 = 30.1; p < 0.0001, inter- 
mediate and high resource concentrations in the same 
homogenous subset, p = 0.11). Furthermore, predator pop- 
ulation density declined over time in all resource concen- 
trations (time: F L9 = 256; p < 0.0001). 

Species proportions in the prey community 

The effects of predation on prey community compos- 
ition differed along the resource concentration axis. The 
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Figure 2 The effects of experimental treatments (resource enrichment and predation) on community composition and population 
densities over time. Vertical bars indicate species proportions in the prey community over time. Horizontal lines represent the total number of 
individuals in the prey community (circles, total number of prey individuals in the prey community) and T. thermophila population density 
(triangles) (mean ±S.E). 



proportion of N. capsulatum was higher than the two 
other species in the low resource concentrations, both 
with and without predators (Figure 2; predators absent: 
F 2 .9 = 58.8; p < 0.0001, S.m and B.c in the same ho- 
mogenous subset, p = 0.93; predators present: F 2 .9 = 35.5; 
p < 0.0001, S.m and B.c in the same homogenous subset, 
p = 0.55). In intermediate resource concentrations, the 
proportions of all species were different from each other; 
the proportion of N. capsulatum was the highest in 
treatments without predators, B. cereus was the lowest, and 
the proportion of B. cereus in treatments with predators 



was the highest and 5. marcescens the lowest (Figure 2; 
predators absent: F 2 .9 = 25.0; p < 0.0001, predators present: 
F 2 .9 = 63.8; p < 0.0001, in both, all species in different 
homogenous subsets). The proportion of N. capsulatum 
was higher than the other two species in high resource con- 
centrations without predators (Figure 2; F 29 = 15.2; p = 
0.001, S.m and B.c in the same homogenous subset, p = 
0.65). The proportion of B. cereus was the highest and 5. 
marcescens was the lowest in high resource concentrations 
with predators (Figure 2; F 2 9 = 116; p < 0.0001, all species in 
different homogenous subsets). 
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Discussion 

Competition for shared resources and predation are usu- 
ally among the most important factors driving prey com- 
munity dynamics [20]. To investigate the effects of these 
two factors on community dynamics after a resource pulse, 
we conducted a microbial microcosm experiment where 
we manipulated the intensity of resource competition with 
resource pulse concentration as well as the presence of 
predation. In order to predict and understand the outcome 
of community dynamics, we measured growth and defense 
related traits of each prey species prior to the long-term 
experiment. Based on these measurements, we were able 
to categorize N. capsulatum as a good resource competi- 
tor, B. cereus being more resistant against grazing and S. 
marcescens as an opportunist specialized in high resource 
environments (Figure 1). Results from the community ex- 
periment support this view: without predators in all re- 
source concentrations and also with predators in low 
resource concentration, N. capsulatum was the dominant 
species (Figure 2). However, predation hindered the domin- 
ance of N. capsulatum; in intermediate and high resource 
concentrations with predation, the inferior competitor but 
more grazing resistant B. cereus was the dominant species 
(Figure 2). See Table 1 for a summary of the results. 

Temporal changes in resource availability, such as re- 
source pulses, are known to promote species coexistence 
[4]. This enhanced diversity is often explained by interspe- 
cific differences in the ability to maintain positive growth 
under low resource conditions and rapid growth when re- 
sources are abundant [21-24]. The role of fluctuating re- 
sources on species co-existence and community dynamics 
was tested in a previous study [17] using the same microbial 
community as here. The main finding in that study was that 
S. marcescens was the dominant prey species across all 
treatments [17], a qualitatively different result compared to 
the present study. However, in the Hiltunen et al. 2012 
study, the inflow of fresh media, having the same concen- 
tration as "high" in the current study, was a frequent event 
(on average -80% of the media was replaced daily), enabling 
S. marcescens to dominant the community. A likely reason 
why S. marcescens was able to dominate in the previous 
study was its high growth rate in the high resource concen- 
tration of the fresh culture media (Figure 1A: high resource 
concentration). In the current study, S. marcescens was 
not the dominant species in any of the treatments. S. 
marcescens seem to be a high resource specialist and in 
our current experiment, even in the initially high re- 
source concentration, the availability of resources de- 
clined rapidly since there was no further input of new 
resources and the consumers depleted the existing re- 
sources. Thus, high resource specialists will experience 
a fitness benefit only for a short period of time and 
after a while species specialized in the use of low re- 
source concentrations, such as N. capsulatum in our 



experiment, will dominate the communities. However, 
we found that in the intermediate and high resource con- 
centrations without predators S. marcescens became the 
dominant species by the last time point and also increased 
its proportion substantially in the low resource environ- 
ment during the last three time points. This is a result that 
cannot be explained with the competitive ability measure- 
ments. One explanation for this change in dominance after 
21 days is the change in the resource quality. Pekkonen 
et al. (2012) studied resource mediated indirect facilitative 
and inhibitory interactions between S. marcescens and N. 
capsulatum, competing for the same plant detritus resource 
as used in our experiment. They found, among other 
things, that S. marcescens benefited from the presence of 
consumed growth media in contrast to N. capsulatum. To 
investigate similar questions, Lawrence et al. (2012) [25] 
used a nuclear magnetic resonance (NMR) spectroscopy to 
study the qualitative changes in a complex, organic plant 
based culture media. This approach could also have been 
useful in our study for investigating qualitative changes in 
our culture media and explain why S. marcescens was able 
to increase its proportion at the end of the experiment. 
However, based on data that we have now, we can conclude 
that in our resource pulse environment, temporal changes 
in the resource availability might not only have been quan- 
titative but also qualitative. Even though we do not have 
direct evidence, this type of complex facilitative interaction 
could explain why the inferior competitor in low concentra- 
tions of fresh growth media could dominate at the end of 
the experiment when concentrations of the original re- 
sources were bound to be extremely low. However, this also 
makes predicting the competitive outcome based solely on 
growth rate measurements more challenging. 

Predation is generally predicted to have a positive ef- 
fect on coexistence among prey species when predators 
prevent the exclusion of more resistant but less competi- 
tive prey types [6,7,16,26]. The community data presented 
here is a good example of this scenario. Without preda- 
tors, N. capsulatum was the dominant species; however, 
with predators in intermediate and high resource concen- 
trations, the good resource competitor N. capsulatum was 
no longer able to dominate the prey community. B. cereus, 
an inferior resource competitor but an inherently more 
grazing resistant species, dominated the prey community 
in these conditions. 

Interestingly, we found an interactive effect between 
resource concentration and predation treatments so that in 
low resource environments N. capsulatum was the domin- 
ant prey species instead of the more grazing resistant B. 
cereus. Holt et al. (1994) were the first to use a theoretical 
approach to investigate the interplay of exploitative and ap- 
parent competition and described the shift from dominance 
of exploitative to apparent competition along the enrich- 
ment gradient. Our experimental findings described above 
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are consistent with the predictions by Holt et al. We ob- 
served, as the theory predicts, a shift in dominance from 
the superior resource competitor (N. capsulatum) to the 
grazing-resistant prey (B. cereus) along an enrichment gra- 
dient when the prey community was exposed to predation. 
Changes in the relative importance of competition and pre- 
dation related traits along an enrichment axis might explain 
this finding in our experiments. When resources are ex- 
tremely limited, traits related to competitive ability are dis- 
proportionally important and the importance of these traits 
may override any benefits that the higher grazing resistance 
can provide. Also, the fact that predator population density 
and, therefore, predation pressure was lower in the low 
resource environment could have contributed to making 
competitive ability a more important factor in determining 
the dominant prey species. 

Conclusions 

We found that in our experimental bacterial communities, 
either resource availability or ciliate predation determined 
the dominant species. This is in line with theoretical pre- 
dictions [1,2,6,7,21,23,26,27]. In most cases, we were able 
to predict the identity of the dominant species based on 
trait characteristics measured prior to the community ex- 
periment in single species short-term assays. We also 
found that resource concentration and predation treat- 
ments strongly interacted so that in low resource environ- 
ments competitive ability was the main factor determining 
the community composition. 

Methods 

Study species and culture media 

Our experimental prey community consisted of three het- 
erotrophic bacterial species: Bacillus cereus (ATCC 14579), 
Serratia marcescens (ATCC 13880) and Novosphingobium 
capsulatum (ATCC 14666), all feeding on shared plant 
detritus medium. All strains were obtained from the 
American Type Culture Collection [28]. The criteria for 
selecting these three bacterial species included positive 
growth on both the nutrient broth agar and liquid prey 
culture medium and that the species were distinguishable 
based on their colony color and morphology. The source 
of organic carbon in our liquid prey culture medium was 
filtered plant detritus (hay extract, Ward's natural science, 
Rochester, NY), for a detailed recipe and methods for cul- 
ture media preparation see [17,18]. Three final concentra- 
tions of plant detritus used throughout experiment were 
0.215, 1.075 and 2.15 mg/1 and they are referred hereafter 
as low, intermediate and high resource concentrations. 
When the experiments were initiated, the bacterial strains 
(stored in 50% glycerol at -70°C) were first thawed and 
then cultivated on nutrient broth agar (10 g nutrient broth, 
2.5 g yeast extract and 15 g agar (all Sigma-Aldrich) in 1 
liter of dH 2 0. 



As a predator, we used a ciliated protozoa, Tetrahymena 
thermophila (ATCC 30008), which is an asexual strain 
consisting of only a single mating type and has been widely 
used in experimental microbial ecology [8,9,17,29-31] T. 
thermophila was cultured in a proteose peptone-yeast 
medium (10 g nutrient broth, 2.5 g yeast extract (Sigma- 
Aldrich) in 1 liter of dH 2 0) in controlled laboratory condi- 
tions prior to the experiment. All experiments were carried 
out at 25 ± 1°C. 

Estimating competitive ability and grazing resistance 

We conducted a short-term growth and feeding experi- 
ment to estimate the competitive ability and grazing resist- 
ance of each prey species. As an estimate of competitive 
ability, we used population growth rate and grazing resist- 
ance estimated as biomass reduction by predators. Mea- 
surements were conducted in each of the three resource 
concentrations (low, intermediate and high). 

When measuring the population growth rates, 350 ul of 
fresh culture medium was inoculated with a small amount 
of bacteria (10 ul containing ~10 A 6 individuals). Then, the 
population was allowed to grow for 96 hours and the initial 
period when resources were not limiting was used to esti- 
mate the maximum growth rate (see below). When mea- 
suring grazing resistance, we used the same bacterial 
cultures that were allowed to reach carrying capacity (96 h) 
before adding 10 ul of centrifuged and washed predator 
stock in pH 7.5 phosphate buffer containing -100 predator 
individuals. We monitored bacterial biomass (optical dens- 
ity) for another 96 hours which allowed predators to reach 
a high density and prey population was simultaneously 
grazed down. This reduction of prey biomass (percent 
grazed biomass from initial biomass prior to the addition 
of the predators) was then used as an estimate of the graz- 
ing resistance. All treatments were replicated nine times. 

Population growth and predator-induced decline was 
measured with Bioscreen C spectrophotometer (Growth 
Curves AB Ltd, Finland) where optical density of each 
well was measured at 480-580 nm wavelengths at five 
minute intervals. Population growth rate was calculated 
as the slope of the linear regression of natural logarithms 
of population biomass versus time when the population 
grew at its maximal rate. This methodology in bacterial 
trait measurements has been used successfully in previ- 
ous studies with similar systems (see e.g. [9,17,18]). 

Microcosm community experiment 

To study community dynamics after a resource pulse, we 
conducted a factorial microcosm experiment in batch cul- 
tures where bacterial communities were cultured in three 
resource levels (low, intermediate and high), with and with- 
out predators. Each treatment was replicated four times. 

The microcosms were contained in 250 ml polycar- 
bonate Erlenmeyer flasks (Corning) and experimental 
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communities were assembled as follows: each bacterial spe- 
cies was grown separately for 96 hours in liquid medium 
and then mixed in even proportions 10 ml each and 750 ul 
of centrifuged predator containing ~ 7.5 x 10 A 4 individuals 
was added to the predation treatments. At the beginning of 
the experiment, the total volume of liquid in each micro- 
cosm was 30.75 ml. The experiment was carried out for 21 
days and sampling was conducted every three days by re- 
moving a 1 ml sample from each microcosm under sterile 
conditions. 

A Lugol's solution was used to fix 200 ul of the 1 ml 
sample and the predator population densities were de- 
termined automatically from the digitized images using 
an image recognition script (for detailed methods see [29]). 
Two predator samples, days 15 and 18, were misplaced 
and treated as missing data. The prey species population 
sizes were determined with serial dilution plating and by 
counting the colony forming units (CFU's) cultured on 
petri dishes containing nutrient broth agar. 

Statistical analysis 

The effect of resource concentration and predation on 
population growth rates and grazing resistance were ana- 
lyzed with a two-way ANOVA. Post-hoc Tukey compari- 
sons were performed to test pairwise interactions. In the 
community experiment, the effect of treatments on popula- 
tion densities and species proportions were analyzed with 
repeated measurements ANOVA (RM ANOVA). The re- 
quirements of the RMANOVA only allowed the analysis of 
the first five samples of the predator population data due to 
the loss of samples from days 15 and 18. Population growth 
rates (h' 1 ) were calculated as r = ln(N t+ i/N t )/t where N t is 
population size at time, t. All analyses were performed with 
PASW statistics (SPSS Inc. Chicago IL, v. 20.0) software. 
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